First passage times and distances along critical curves 
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We propose a model for anomalous transport in inhomogeneous environments, such as fractured 
rocks, in which particles move only along pre-existing self-similar curves (cracks). The stochastic 
Loewner equation is used to efficiently generate such curves with tunable fractal dimension d/. 
We numerically compute the probability of first passage (in length or time) from one point on the 
edge of the semi-infinite plane to any point on the semi-circle of radius R. The scaled probability 
distributions have a variance which increases with df, a non-monotonic skewness, and tails that 
decay faster than a simple exponential. The latter is in sharp contrast to predictions based on 
fractional dynamics and provides an experimental signature for our model. 



INTRODUCTION 

Scale invariant curves, such as coastlines and fracture 
fronts, abound in nature [ij. The fractal dimension df 
relates the scaling of the length £ along the curve between 
two points to their actual separation R in Euclidean space 
hy i (X R'^f . However, for stochastic shapes, the ac- 
tual length will fluctuate, and knowledge of its full dis- 
tribution greatly augments the information from scaling. 
We numerically investigate some related distributions for 
scale invariant (critical) curves in two dimensions. 

While we measure the distribution of lengths along 
fractal curves, an underlying motivation is to understand 
the distribution of times for traversing such curves. In 
particular, we are interested in exploring flow through 
porous media, such as diffusion of chemical tracers in- 
jected into a water-saturated heterogeneous permeable 
sediment or fractured hard rock formation. Experiments 
suggest that flow proceeds along a network of pre-existing 
fractal paths [2] , as further corroborated by simulations 
modelling streams as percolation clusters Q. In such 
cases, flnding the distribution of first passage times be- 
tween points of known Euclidean distance requires knowl- 
edge of the dynamics of the process, as well as the dis- 
tribution of traversed lengths. If we assume that the 
motion along the path is independent of its shape (as in 
a simple damage spreading process, or in curvature in- 
dependent diffusion), we obtain the distribution of first 
passage times by simple manipulation of the distribu- 
tion of lengths. At the minimum, results based on such 
an assumption help to rule out simple hypotheses about 
the dynamics, and/or structures of the fractal shapes, in 
physical models. 

Important examples of scale invariant curves in two 
dimensions are self-avoiding walks (SAWs) and percola- 
tion fronts (PFs); additional instances may be obtained 
from interacting systems at a critical point These 
examples suggest that generating and studying stochas- 
tic scale-invariant curves is a computationally hard pro- 
cess. Indeed, the past decades have witnessed much effort 



and progress in efficient numerical algorithms for gener- 
ating SAWs and critical systems @ . An important recent 
development concerns the Stochastic Loewner Evolution 
(SLE), which provides a way of generating such curves 
through mappings (in the complex plane) of a simple 
Markovian random walk 0,0, 0- While the main inter- 
est in SLE has been to extract analytic information about 
critical systems, here we use it as an efficient means to 
generate curves with tunable fractal dimension. 



SLE, CRITICAL CURVES AND THEIR LENGTHS 

Specifically, we study non-self-intersecting curves in 
two-dimensions that begin at the origin and remain in the 
upper half-plane. Such curves can be parameterized by a 
complex function 7(i) of a real argument t. Loewner has 
shown that for any such curve described by j(t'), with 
< t' < t, there exists a conformal mapping gt{z) that 
maps the upper half-plane surrounding the curve into the 
complete half-plane, while the curve itself is mapped onto 
the real axis .10,] . In particular, for the tip of the curve 
Qtilit)) = C(i)j where ((t) is a real function. According 
to Loewner, the mapping gt satisfies the equation 



dgtjz) 
dt 



gt{z)-at)' 



(1) 



This equation has a natural parameterization in terms 
of the so-called Loewner time t, which has dimensions 
of an area and whose definition is based on a standard 
form of 01 ■ The knowledge of the real function C(t) 
enables a complete reconstruction of the curve j{t): this 
is accomplished by integration of Eq. ([T]) backwards in 
time, starting from the "final condition" (?t(7(i)) = C{t) 
and using gt=o{z) = z. 



It has been shown 0, [1, S 



11 , 1^ that different sets of 
scale invariant curves with fractal dimensions 1 < rfy < 
2 map onto standard one-dimensional Brownian motion 
functions C(0 — \/^B(t), with zero mean and with a 
variance k related to the fractal dimension by d/ = 1 + 
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FIG. 1: A discretized implementation of SLEg/3, which in the 
scaling limit is conjectured to be a SAW in the upper half- 
plane. The curve is generated by integrating Eq. ((T)) backward 
in time with step size dt = 10~^. The trace j{t) is stopped 
upon its first encounter with the semicircle of radius R — 1. 
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k/8 [13|. This process is denoted as SLE^ and can be 
used to conveniently generate long scale-invariant curves. 
For particular values of k the traces "/(t) coincide with 
the continuum limit of well known curves, such as SAW 
or PF, defined on critical lattice models [8!]. An example 
curve generated by implementing SLEg/a (corresponding 
to a SAW) with discrete time steps is shown in Fig. [TJ 

For SLEk traces such as in Fig. [U an important is- 
sue is how to measure length between two points on the 
curve. In fact, continuous-time Brownian motion gener- 
ates SLEk paths which contain infinitely minute details, 
thus having an infinite length. This problem has been 
addressed by Kennedy [3, and like him we intro- 
duce a specific procedure for calculating the lengths of 
curves generated by time-step discretization. Loewner 
time is discretized into m — t/ dt equal time steps. While 
the discrete time points idt are equally spaced, Loewner 
equation produces a discrete set of positions z(idt), where 
the two-dimensional distances between sequential points 
Ai = \z(idt) — z{{i — 1)^^)1 are extremely non-uniform. 
For the solution of the physical problem we need to know 
the length of the walk between such points. In the lat- 
tice version of the walk, we expect the number of steps 
between two points to be proportional to A'*^ . One can 
use this relation to define the length £ of the walk as 



FIG. 2: Comparison between the rescaled distributions of 
lengths pdf{(./R^f) for SLEk and for the corresponding lat- 
tice models. The solid line with a narrow peak corresponds to 
the probability distribution of SLE with k = 8/3 and is com- 
pared with that of SAWs, plotted in crosses. The solid line 
with a wide peak corresponds to the probability distribution 
of SLE with K = 6 and is compared with that of PFs, plotted 
in dots. All distributions have been further rescaled so that 
their average is equal to one. 



Eq. ^ converges to a finite limit. Since we use the defi- 
nition in the context of first passage problems, applied to 
a stochastic variable, £, which is diff'erent from those con- 
sidered in 14 , 1^ , it is essential to verify the correctness 
of the procedure. We do so by comparing distributions 



obtained from lattice implementation of SAWs and PFs 
with those obtained from SLE„ at corresponding values 
of K. 



DISTRIBUTIONS OF FIRST PASSAGE 
LENGTHS 



1=1 



Relying upon Eq. we have computed the distribu- 
tion of lengths for SLE^ with k = 8/3 and k = 6 for 
curves starting from the origin of the upper half-plane 
and arrested upon touching for the first time a semicircle 
of radius R. For these values of k, SLE^ is conjectured 
and shown to represent the scaling limit of SAWs on a 
square lattice _and PFs on a triangular lattice, respec- 
factor c with dimensions of tively 17, 



£ EE cJ2 l^i^dt) - z{{i - l)dt)ff, 



(2) 



with a proportionality 

Consistency and suitability of this definition for the 
standard relations between the number of steps and the 
end-to-end distance of the curve and its proper behavior 
in the m ^ 00 limit have been demonstrated numerically 
by Kennedy IJ, [15|. As dt 0, more points are sam- 
pled, their separations become smaller, but the sum in 



18l |. We have computed via Monte Carlo 
simulation the length distributions for these lattice mod- 
els (which have a natural definition of length), adopting 
the same boundary conditions. In Fig. [2] we compare 
the rescaled pdf{£/R'^f) for the two sets of distributions: 
the numerical value of c in Eq. [2] for SLE^ is fixed by 
imposing that the average of the distribution is one. 
The agreement between SLE„ and the corresponding 
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FIG. 3: Rescaled distributions of lengths for SLE^ with k — 
0.01, 0.1, 0.2, 0.5, 1, 1.5, 2, 2.5, 3, 3.5 and 4. The distributions 
have been further rescaled so that their average is one. As 
K (df) increases, the peak of the distribution moves to the 
left while its width becomes larger. Dashed lines represent 
the rescaled distributions for PFs (corresponding to k = 6) 
for R = 100 and R = 500. It should be noted that the two 
distributions are almost indistinguishable, and clearly show a 
faster than exponential decay. The inset shows the skewness, 
X3, of the length distributions as a function of k. 

lattice models is quite good, thus supporting our choice 
of Eq. Q as a proper means of measuring lengths. The 
quantity which is calculated as £ is therefore the true 
length (in the usual sense) of a polymer chain for k — 8/3, 
of a percolation front for k — 6, or in general of a scale 
invariant self-avoiding curve. The somewhat unintuitive 
form of Eq. ([2|) is a consequence of Loewner evolution, 
which, when used as a finite difference equation with 
equally spaced time instances U, produces a non- uniform 
distribution of Euclidean distances z(ti + dt) — z(ti) [Tst - 

However, we should stress that there are also notable 
differences between the two sets of curves: one possible 
source of discrepancies comes from finite-size (lattice) ef- 
fects, since SLE^ corresponds to a continuum limit in 
which the number of steps on a lattice goes to infinity, 
while the lattice spacing goes to zero. A second source is 
intrinsic to SLE^ for k > 4. At these values the traces can 
touch themselves, the numerical integration of Eq. ([T]) be- 
comes problematic fis'l, and convergence to the asymp- 
totic limit is slow. 

Having gained confidence in the validity of our numer- 
ical measure of length along critical curves, we use SLE^ 
to generate curves for a range oi df for which there is 
no simple lattice analog. It is worth emphasizing that 
the appearance of SLE^ enables dealing with problems 
which were previously untreatable: although SAWs and 
PFs can also be calculated directly, most other scale- 
invariant curves with a generic df cannot be produced 
by any other method. The resulting scaled probability 



distributions of the lengths of these curves, pdf{i/R'^^), 
are plotted in Fig. [3] Each distribution is obtained from 
10^ SLE„ trajectories starting from the origin and ar- 
rested upon first touching a semicircle with radius R. 
The integration step for all curves is dt = 5 x 10"'*. The 
distributions depend only on the scaled parameter £/R''-f , 
and we have further rescaled both axes such that the 
mean value is at unity. The distributions are asymmetric 
and become broader as d/ is increased. The variance ini- 
tially grows quadratically in k, then levels off to a slower 
growth. Interestingly, the skewness (a measure of asym- 
metry) is not monotonic in d/ and has a maximum at 
K « 1, as depicted in the inset of Fig. [3] The entire 
shape of the distribution (and by extension the variance 
and the skewness) are universal functions of df. 

For comparison with other systems, an important ques- 
tion concerns how the distribution behaves at large ar- 
guments. Detailed numerical simulations of PFs allow 
us to conclude that the tails of the distributions decay 
faster than exponentially. We generated PFs for the rel- 
atively large sizes of i? = 100 and R — 500 (10^ and 
10^ trajectories, respectively). As indicated in Fig. [31 a 
non-linear behavior at large arguments is present in this 
semi-log plot, irrespective of R (implying that finite-size 
effects can be neglected). Intuitively, we expect this fea- 
ture to arise from the self- avoiding nature of SLE^: for 
a Markovian path allowed to intersect itself, the prob- 
ability of adding a new segment to an existing trace in 
the domain is independent of its length. Such paths are 
characterized by a distribution that decays as a simple 
exponential. For SLE^, by contrast, the available space 
within the domain decreases with increasing t {£) and the 
escape probability increases. 

APPLICATION TO TRANSPORT IN POROUS 
MEDIA 

Our results lead to a practical application in the con- 
text of anomalous transport: as a specific example, we 
focus on the paradigmatic case of flow in rocks, in which 
motion is constrained to proceed along preexisting (frac- 
tal) curves. Non-diffusive dynamics is quite common 
in inhomogeneous environments such as these, and fre- 
quently modeled by a fractional Fokker-Planck equation 
(FFPE) [ll, [lOl. In the FFPE approach, the medium 
where particles move is supposed to be globally homo- 
geneous, and the effects of heterogeneities are taken into 
account by assuming that particles are locally trapped 
and thereby experience extremely long (power-law dis- 
tributed) waiting times. This is then mathematically 
modelled through fractional derivatives in time. How- 
ever, the behavior resulting from specific physical real- 
izations (such as motion along scale invariant curves) can 
be quite different from that arising from ideal fractional 
dynamics. To elucidate these differences, we focus on 
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FIG. 4: FPTD along critical curves with df = 1.5 (a); con- 
trasted with the results for normal diffusion (b); and FFFE 
(c). Left: diffusion-dominated transport (Pe=0). Right: 
drift- diffusion (Pe=10). Time is in dimensionless units f = 
tD/E? . Note that for (c) the mean first passage time is infi- 
nite. 



the first passage time distribution (FPTD), namely the 
probability that a walker starting at the origin in Fig. [1] 
reaches the circle at radius R at time t for the first time. 

Within the FFPE approach, the trapping events lead 
to the anomalous transport scaling r oc i?^, with z > 2. 
Analysis of the equation then indicates that for a domain 
of finite size the FPTD has an asymptotic power-law 
decay as r~^~^/^ for times much longer than the typi- 
cal time scale of the process (the specific nature of the 
boundaries does not affect this scaling) 23| 2^- FFPE 
also allows incorporating an additional external drift ve- 
locity u, and in this case the FTPD (after an initial tran- 
sient whose duration is proportional to the ratio R/v) 



decays as t^^^^/^ , with z' > 1 [2l|[2^. In both cases, 
the long retention times are such that the mean time to 
reach the outer boundary is infinite 20]. 

In our model of flow along scale-invariant streams, the 
time it takes to first reach a Euclidean distance R de- 
pends on two elements. One is the dynamics by which 
motion proceeds along these tracks, and the second is the 
distribution of lengths £ of appropriate tracks for a given 
Euclidean distance R. The latter is the quantity that 
we calculated in Fig. O the former can depend on many 
features of the curve (such as curvature, vicinity of other 
portions) as well as on its length. To make the problem 
tractable, we make the simple assumption that the prob- 
ability that the walker reaches its destination (for the first 
time) at a time t only depends on the length £, and is 
given by some function p(r|£). The overall FPTD is then 
obtained by the integral pdf{T\R) = J dip{T\£)pdf{i\R)- 

In a damage spreading process (such as burning fire) 
in which the walker proceeds at a constant velocity, 
p{t\£) = 5{vT — £). More generally, for a drift-diffusion 
process we must solve the eigenvalue problem for first 
passage over a segment of length £, which depends on the 
Peclet number Pe= vi/D and on the initial and bound- 
ary conditions on the interval ^22i] . To be consistent with 



the macroscopic boundary conditions on the semicircle, 
we impose that the particles are injected at the origin 
and can escape from the track only upon touching the 
other extreme of the curve at the outer boundary of the 
traversed medium, where they are absorbed. Mathemati- 
cally, this is implemented by assigning a reflecting bound- 
ary condition at the origin (any particles moving back to 
the injection point are sent back without loss), and an 
absorbing boundary on the end point at £. Under these 
assumptions, the first passage time distribution reads as 
follows: 



p{T\£)^2e^J2- 



1 



C 



Xj sin(Aj), 



(3) 



where Xj are the roots of tan(A) = — A/C, and ( = Pe/2 

In Fig. m we plot the resulting FPTD, pdf{T\R), for 
particles diffusing along SLE„ (obtained as a numerical 
convolution). The left figure corresponds to pure dif- 
fusion (Pe — 0) along the curve, which appears sub- 
diffusive in Euclidean space, since t (x £'^ (x R^'^f and 
thus z = 2df > 2. The results are contrasted in the 
same figure with normal diffusion (exponentially decay- 
ing FPTD), and with that of FFPE (power-law decay). 
FPTDs for finite values of Pe are plotted on the right 
side, which appear as super-diffusion in real space, as 
T (X £ (X R'^f , whence z' = df < 2. The corresponding 
curves for normal diffusion and FFPE are also shown. 
As discussed before, when motion along the curve pro- 
ceeds at uniform velocity (no diffusion and Pe — > c») 
the tails of the distribution decay faster than a simple 
exponential. Thus, in all cases the presence of a fairly 
wide distribution of path lengths broadens the FPTD, 
but the tails of the distribution fall off more rapidly than 
predicted by fractional dynamics: this allows for a very 
strong discrimination between other self-intersecting and 
non-self-intersecting trajectories at microscopic level. (In 
particular the MFPT is finite.) This difference is not de- 
tectable by just looking at the macroscopic anomalous 
transport scaling exponent. Moreover, since FPTD is 
usually an easily experimentally accessible quantity, this 
observation may provide a way to support or reject the 
hypotheses about the microscopic dynamics which lie at 
the basis of these models. 



CONCLUSIONS 

In this Letter we have proposed a model for anoma- 
lous transport in inhomogeneous environments, which 
is suitable to capture, e.g., the dynamics of flow 
streams through fractured rocks. In particular, the two- 
dimensional probability density of the particles which can 
be obtained from our model qualitatively agrees with the 
typical distribution of a chemical plume flowing through a 
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two-dimensional fracture network (as shown in Fig. 4(b) 
ofRef. [H). 

Unlike fractional dynamics, the FPTD docs not fall off 
as a power law and has a finite mean. A necessary ele- 
ment of the model is the distribution of lengths of fractal 
curves, which we obtain by numerical studies of lattice 
models and SLE: to the best of our knowledge, the first 
passage distributions obtained here have not been de- 
rived analytically. We hope that analytical approaches 
to SLE can provide exact insights on these distributions, 
and hence enable quantitative comparisons with experi- 
ments. 
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